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Metal-Fluorocarbon Pyrolants. XV: Combustion of two 
Ytterbium—Halocarbon Formulations 


Ernst-Christian Koch,**” Volker Weiser,‘ Evelin Roth,‘ Sebastian Knapp,*‘ 


Joost van Lingen‘ and Jacco Moorhoff* 
“NATO Munitions Safety Information Analysis Center (MSIAC), Boulevard Leopold II, B-1110 Bruxelles, 
Belgium 
> Technische Universitat Kaiserslautern, Erwin-Schrédinger-Str. 52, D-67663 Kaiserslautern, Germany 
° Fraunhofer Institut fiir Chemische Technologie, Joseph-von-Fraunhofer-Strasse 7, D-76327 Pfinztal, Germany 
4 TNO Technical Sciences, Lange Kleiweg 137, 2288 GJ Rijswijk, Netherlands 


Abstract: Binary pyrolants based on ytterbium and either polytetrafluoroethylene or hexachloroethane 
yield a broad stable combustion regime with intense luminous flames. Combustion temperatures have been 
determined spectroscopically for both pyrolant types. The temperature correlates with the burn rate of the 
pyrolant. The UV-Vis spectra of the fuel rich pyrolants have been measured and were converted into x,y- 
color values in the 1931 CIE color diagram. The color obtained with both pyrolant types produce shorter 
dominant wavelength emission and higher saturation than common BaCl BO) emitters. 


Keywords: Ytterbium, Polytetrafluoroethylene, PTFE, Pyrolant, Burn rate, Combustion Temperature 


Introduction 


Some time ago it was reasoned by Sturman that 
the rare earth metals and their compounds should 
in principle function as flame color agents.’ The 
reason for his assumption was the observation that 
the individual elements give distinct emissions 
when introduced as aqueous salt solutions into 
inductively coupled plasma. Earlier Dumont 
had found that lanthanum acetylacetonate when 
used as fuel in an ammonium perchlorate based 
formulation yields a white flame.” Mischmetall and 
cerium acetylacetonate when used in perchlorate 
based formulations both yield a faint orange flame.* 
Another important prediction made by Sturman 
was that either the rare earth metals or their alloys 
should also yield sparks with colours other than 
those typically found on the black-body curve. 
That is, distinct red and green emissions should 
be accessible at best through use of some of these 
elements. However apart from the above work 
the potential of rare earth metals in pyrotechnic 
formulations has never been investigated. 
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Recently we have shown that the lanthanide metals 
samarium (Sm), europium (Eu), thulium (Tm) 
and ytterbium (Yb) can burn in the vapor phase.* 
This is due to compliance of these metals with 
Glassman’s criterion which requires vaporization 
temperatures of the metals to be lower than the 
dissociation temperatures of the corresponding 
oxides.*° Powdered ytterbium undergoes vapour 
phase combustion in oxygen and with halogen 
based oxidizers.° A narrow range of ternary 
mixtures of Yb with polytetrafluoroethylene and 
Viton® yields extended luminous flames which 
radiate intensely in the infrared spectral range and 
outperform standard infrared decoy flare payloads 
based on Magnesium/PTFE/Viton® (MTV)” in a 
certain spectral band.* In view of the unexpected 
results of ytterbium based formulations we decided 
to investigate the combustion behaviour of binary 
mixtures of ytterbium with polytetrafluoroethylene 
and hexachloroethane, the latter being the preferred 
oxidizer in obscurant formulations.’ 
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Experimental details 
Warning 


Metal-fluorocarbon  pyrolants are explosive 
materials inherently sensitive to electrostatic 
discharge, heat, flame and impact. These materials 
when alight yield intense radiation sufficient to 
cause severe skin burns. In the unconsolidated 
state the material can undergo deflagration and 
even shock up to low order detonation in larger 
quantities (>1 kg).”° 


All preparations must be conducted in accordance 
with appropriate national safety regulations. In 
particular the personnel involved should wear 
flame-resistant personal protection equipment 
such as overalls and balaclava made from 
NOMEX® -III or aluminized PBI® with a facial 
heat shield. The personnel must be grounded by 
conductive protective shoes and a wristband, and 
all equipment and tooling must be grounded as 
well.” 


Ytterbium ingot and powder (99.9% rare 
earth oxide (REO)), volume median diameter 
D{v, 0.5] = 3.19 um, polytetrafluoroethylene, 
D{y, 0.5] = 7.5 wm and hexachloroethane were 
all purchased from ALFA AESAR and used 
without further treatment. The ytterbium ingot 
was rasped down with a rasp to yield a coarse 
powder of ~150 um mean diameter. 5 g batches of 
the pyrolants were prepared by dry mixing of the 
components in a 25 ml spherical mixing container 
made of conductive polypropylene for 30 min at 
120 rpm in a tumbling blender. All pellets were 
pressed at 10 MPa pressure applied for 10 s ina 
5.6 mm cylindrical die under safety precautions. 
Typical sample weight ranged between 1800 
and 2000 mg. Pellets were firmly wrapped in 
transparent adhesive Tesa®-film to protect the 
lateral surface from flame spreading but to keep 
the burn front visible. The cylinders were glued 
with cyanoacrylate glue on small ceramic squares 
(20 x 20 x 4 mm). The pellets were transferred 
into a fireproof fumehood with constant venting 
and ignited with a non-luminous propane torch 
flame. 


For the temperature determination Near infrared- 
(NIR) spectra were recorded with a PGS-NIR- 
Spektrometer 2.2 (Carl Zeiss AG, Germany) in the 
wavelength range 2 = 1.0—2.1 um and calibrated 
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with a technical blackbody heated to 2000 K. Using 
proprietary ICT-BaM-Code calculated spectra 
were fitted to the experimental spectra with a least 
squares fit.'° The UV-Vis spectra were recorded 
with a spectrograph Shamrock A-SR-500i-B2 with 
Andor DU920P-UV-DD-detector (both purchased 
from Andor Inc. USA). A diffraction grating with 
150 lines/mm was used. 


The combustion process was recorded with a 
Panasonic HDC-HS20 with the following settings; 
aperture: 16; exposure time: 1/8000 s. The 
accuracy of the burn rate determination is limited 
by the frame rate of the video camera which is 
28 frames/second. 


Results and discussion 


Pressing densities (p,,,) greater than 90% 
Theoretical Maximum Density (TMD) were 
achieved for Yb/PTFE-compositions containing 
up to 82 wt% Yb. At higher ytterbium contents the 
mixtures do not consolidate in the same manner 
and yield densities of around 87% TMD (Fig. 1). 
Stable combustion of Yb/PTFE is observed from 
62-95 wt% Yb. In this range the combustion rate 
increases exponentially from about 1 mm s ' to 
26 mm s | at &Yb) ~ 90 wt% and 87% TMD. 
With further increasing Yb content the combustion 
rate decreases again to around ¢(Yb) ~ 95 wt% 
as is depicted in Fig. 2. Below 62% Yb ignition 
of the pyrolant occurs but is followed by rapid 


] 100 % TMD 


_- 90 % TMD 


Density, p [g cm"] 


nr | er er er ee nr ae 
Mass Fraction, &(Yb) [wt-%] 


Figure 1. Experimental density of ytterbium/ 
polytetrafluoroethylene pyrolants and theoretical 
maximum density (TMD). 
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Figure 2. Combustion rate of Yb-based pyrolants 
and combustion temperature with bars indicating 
the scatter of measured data. 


extinction. Above ~95 wt% Yb no cigarette-type 
burn is observed but the reaction proceeds with 
very fast glowing of the complete strand followed 
by delayed scattering of the material, probably due 
to afterburn with atmospheric oxygen. Thus no 
proper propagation rate could be determined. The 
burn rates for binary Yb/PTFE are higher than for 
three ternary Yb/PTFE/Viton pyrolants reported 
recently.* This is mainly due to higher porosity 
(eqn 1) of the present samples (yg = 10-20% as 
opposed to 1—2% for the samples prepared in 
Ref. 8) which facilitates heat transfer via filtrating 
combustion. '! 


_ Pexp 


1 
TMD 2 


g=1 


The spectroscopically determined combustion 
temperature is superimposed for 65, 70 and 80 wt% 
Yb each on the burn rate in Fig. 2. The measured 
mean temperatures (1979, 2000 and 2250 K) 
parallel the combustion rate. A similar behavior 
has been observed with other metal—fluorocarbon 
pyrolants recently.'*'? 


Pressing densities (p,,,) between 80-90% TMD 
were achieved for Yb/HC-compositions (Fig. 3). 
Stable combustion of Yb/PTFE is observed from 
55-85 wt% Yb. In this range the combustion rate 
increases exponentially from about 2 mm s ! to 
17mms ! at &(Yb)=75 wt% and 85% TMD. With 
further increasing Yb content the combustion rate 
decreases again to around €(Yb) ~ 85 wt% as is 
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Figure 3. Experimental density of ytterbium/ 
hexachloroethane pyrolants and theoretical 
maximum density (TMD). 


depicted in Fig. 4. 


Below 55% Yb ignition of the pyrolant 
occurred, but was followed by rapid extinction. 
Above ~85 wt% Yb no cigarette-type burn is 
observed, but the reaction proceeded with very 
fast glowing of the complete strand followed 
by delayed scattering of the material, probably 
due to afterburn with atmospheric oxygen. 
The spectroscopically determined combustion 
temperature is superimposed for 59 and 75 wt% 
Yb each on the burn rate in Fig. 4. The measured 
mean temperatures (1700 and 2100 K) nicely 
parallel the combustion rate. 
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Figure 4. Combustion rate of Yb-based pyrolants 
and combustion temperature with bars indicating 
the scatter of measured data. 
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Fuel rich formulations, €(Yb) = 70 wt% of both 
Yb/PTFE and Yb/HC yield a distinct green 
coloration of the flame that however is partly 
overlapped by lateral orange soot emissions. 
Thus the UV-Vis spectra of both pyrolants yield 
a blackbody continuum superimposed from 
molecular emissions of Yb species Yb, YbO 
YbF and YbCI respectively as shown in Fig. | of 
Ref. 6. Subtraction of the blackbody curve yields 
the principal selective emitters as shown in Fig. 5 
and 6. 


Normalisation of the spectra, separate convolution 


2.05 


Intensity, I, [a.u.] 
oO 
ii 


Wavelength, 4 [nm] 


Figure 5. UV-Vis composite spectrum from YbF, 
YbO and Yb. 
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Fig. 6 UV-Vis composite spectrum from YbCI, 
YbO and Yb. 
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with each of the three tri-stimulus color matching 
functions,'* and integration yields the XY, Y and 
Z-values which were converted to x and y-values as 
shown in Table 1| and in Fig. 7 and Fig. 8. Dominant 
wavelength, A, given by ©, and saturation, », were 
determined graphically. 


Due to very similar electronic transitions the 


Figure 7. Color locus of Yb, YbF; YbO composite 
spectrum,+, dominant wavelength, ©, and white 


point, x. 


Figure 8. Color locus of Yb, YbCl, YbO composite 
spectrum,+, dominant wavelength, ©, and white 
point, x. 
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Table 1. Spectral properties of the composite emitter spectra based on Yb and reference data taken from 


ref. 19 
Yb/PTFE Yb/HC Ba(NO,),/Mg/PVC* B/KNO,* B,C/KNO,* 
Aq (nm) 537.6 541.2 562.3 559.3 561.9 
2 (%) 69.0 73.8 61.5 55.0 52.0 
* plus organic binder for each formulation 
emission spectra of YbCl and YbF'*"’ are very upon ignition with a propane torch yields a loud 


similar and hence close color coordinates are 
obtained. The calculated dominant wavelengths 
are 2 = 537.6 and 541.2 nm, the saturations are 
69% and 73% respectively. 


For comparison Table 1 shows data for a 
conventional green star based on barium nitrate 
(Ba(NO;),)/magnesium/PVC/binder, boron(B)/ 
potassium nitrate (KNO;)/binder and _ boron 
carbide (B,C)/KNO,/binder'* *° and the tested Yb/ 
HC and Yb/PTFE. Yb based formulations produce 
both lower wavelength and more saturated green 
than common formulations based on barium 
nitrate boron or boron carbide. This underlines the 
importance of investigating the rare earth metals as 
color agents in pyrotechnic flames. Investigation 
of organic ytterbium salts as fuels in smoke free 
composition is underway. 


Consolidated pellets of freshly rasped coarse 
ytterbium (150 um) mixed with PTFE (86/14 wt%) 


crackling sound and nice green sparks. A closer 
view actually shows that the actual sparks are 
white but are surrounded by areas with green 
chemiluminescence (Fig. 9). 


Conclusion 


Ytterbium/polytetrafluoroethylene and ytterbium/ 
hexachloroethanepyrolantsyieldstablecombustion 
regimes with intense luminous flames ranging 
from 62—95 wt% and 55-85 wt% Yb respectively. 
The combustion temperatures derived from NIR 
spectra for both pyrolant types correlate with the 
burn rate. The principal emitters in the UV-Vis 
range have been extracted and converted into x,y- 
color values in the 1931 CIE colour diagram. The 
colors obtained with both pyrolant types produce 
shorter dominant wavelength emission and higher 
saturation than common BaCl or BO, type emitters. 
Combustion of pyrolant using coarse ytterbium 
yields intense green sparks. 


Figure 9. Crackling combustion of Yb/PTFE (86/14) ejecting green sparks. 
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Numerical Simulation and Validation of Pyrotechnic Smoke 
Emissions 


Abdelkarim Habib and Christian Lohrer 
Federal Institute for Materials Research and Testing (BAM) - Unter den Eichen 87, 12205 Berlin, Germany 
Email: christian.lohrer@bam.de 


Abstract: This paper aims at evaluating the practicability of an alternative strategy for the classical 
assessment of possible hazards aligned with the use of theatrical pyrotechnics especially at indoor venues, 
where the emission of solid reaction products (aerosols) is an important factor for permission for the use 
of pyrotechnics at stage shows. A CFD model was used to calculate the respective aerosol liberation and 
dispersion in a specific indoor facility. The results from these numerical simulations were compared with 
respective aerosol number concentration measurements carried out before, during, and after the burn- 
off of theatrical pyrotechnic fountains in an examination room, representing a small stage or theatre. 
The simulation results reveal a fair agreement with the experimental results in general, with the trend to 
conservatively over-predict the aerosol concentrations at close range to the pyrotechnic article. The time 
dependent increase and decline of the aerosol concentrations due to the ventilation conditions as observed 
in the preceding experimental studies could be reproduced in the simulations. 


Keywords: Theatrical pyrotechnic articles, aerosols, CFD 


In general, pyrotechnic articles to be placed on 
the European market fall under the scope of the 
Directive 2007/23/EC.' Theatrical pyrotechnic 
articles are categorized in the following 
categories: 


Introduction and background 


The use of pyrotechnics, such as theatrical articles, 
for numerous stage and television productions, 
is globally increasing. New effects and further 
developments widen this field of application 
almost every week. Due to this challenging 
evolution, effective measures are necessary in 
order to assess the possible risks aligned with 
the use of such articles, especially during indoor 
productions, where the audience and the production 7 
staff involved are exposed to potentially harmful 
reaction products of relevant concentrations. 
Enforcement bodies, e.g. those responsible for the 
permission for shows including the use of various 
pyrotechnic articles, are often overwhelmed by 
the large variety of different special effects and 
the impacts of the reaction products associated 
with them. This leads in many cases to a complete 
rejection of the use of pyrotechnics at a specific 


e Tl: pyrotechnic articles for stage use which 
present a low hazard, with a further sub- 
categorization T1 “for outdoor use only”, if 
considered necessary. 


T2: pyrotechnic articles for stage use which 
are intended for use only by persons with 
specialist knowledge. 


The provisions of that directive assure that these 
articles are first introduced to a notified body in 
order to get them appropriately type tested and to 
demonstrate that all essential safety requirements 
of this Directive are fulfilled. Furthermore, the 
articles produced after that EC type-examination 


venue, just to avoid a situation where the enforcers 
could be blamed for not acknowledging all possible 
hazards. 


Article Details 
Manuscript Received:- 28/08/2012 
Publication Date:- 20/10/2012 


Page 10 


must also be in conformity to the initial type. This 
is usually guaranteed by applying sufficient quality 
systems for production and/or end product testing 
(e.g. batch testing), which are subject to regular 


Article No: - 0094 
Final Revisions:-20/10/2012 


Archive Reference:-1506 
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inspections by notified bodies, as well. 


However, the technical requirements for certif- 
ication purposes as given in the applicable standard 
series prEN 16256° cover only performance 
parameters (such as effect dimensions, burning 
rate and sound pressure level etc.), net explosive 
content (NEC) and properties of the articles before, 
during, and after functioning (integrity, stability, 
unintended explosions, burning matter on ground 
etc.). Explicit requirements on gaseous or solid 
reaction products are not included in this standard. 
The categorization into T1 “for outdoor use only” 
may only be based on subjective observations of 
aerosol emissions during the EC type-examination 
tests by the staff, if all other requirements are met. 
Furthermore, how does a parallel use of numerous 
articles, in which all single articles have an indoor 
use permission, change the authorization of the 
use of such articles in a theatrical show? A reliable 
forecast of possible hazards due to the parallel 
firing of diverse articles is almost impossible 
for the pyrotechnicians and the enforcement 
bodies. Confident permission for an envisaged 
use of theatrical pyrotechnic articles can only be 
given if all relevant boundary conditions, such as 
ventilation systems, room geometry or any other 
protective measures of the affected persons are 
considered. 


eles 


Numerical simulations, e.g. computational fluid 
dynamics (CFD), may play a major key role 
in future evaluations of such cases, as they can 
be considered as a possible alternative in this 
application field to estimate the potential hazard 
of the use of specific pyrotechnics under defined 
conditions. CFD is an effective and powerful 
tool, especially when considering the increasing 
development of computational processing power 
in the recent past. In addition, the results from 
CFD calculations might be integrated in quantified 
risk assessments (QRA), as well, in order to get a 
more reliable statement. 


The objective of this work was to investigate the 
practicability of CFD in terms of accuracy of 
results, time and effort, compared to conventional 
estimations of the impacts of solid reaction 
products (aerosols) on a possible indoor use of 
the respective pyrotechnic articles. This has to be 
seen in the context of making the authorization 
process easier for the involved parties: applicants 
(pyrotechnicians or manufacturers) and 
enforcement bodies. 


For model validation purposes, the respective 
CFD calculations of this work, using ANSYS CFX 
V14.0, were compared with the results obtained 
from previous aerosol measurements during the 
burn-off of theatrical pyrotechnic articles at an 


Figure 1. Examination room for the measurements carried out by Dutschke et al. 
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indoor venue by Dutschke et al.* 


Previous experimental setup — 
overview 


In order to get a reliable validation of the applied 
numerical model, a case study was chosen in which 
aerosol emissions during the burn-off of theatrical 
pyrotechnic articles in a venue-like environment 
were investigated. These results were previously 
published by Dutschke et al.* In this work the 
authors measured particle number concentrations 
in the indoor air during and after the burn-off 
of fountains (amongst others) under defined 
ventilation conditions. The experiments were 
performed in a lecture hall with several seating 
rows, depicting a theatre or stage, see Figure 1. 


The room had an overall volume of approximately 
718 m°, offering the following ventilation 
conditions: 


* Outlet—integral: 4400 m*h ' (directly through 


Distance to source Point: 
Point 1: Ax=0m; Ay=3m; Az=1.5m 
Point 2: Ax=6m; Ay=0m; Az=3.1m 
Point 3: Ax=10m; Ay=0m; Az=4.7m 


Figure 2. CAD model of the examination room. 
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openings at the room ceiling), and 


* Inlet — integral: 4200 m? h' (directly through 
circular openings underneath each seat and 
implicitly through the clearance between 
both door frames and door panels due to the 
dominating outlet ventilation conditions). 


The theatrical fountains were burned off centered 
between both laboratory benches, see Figure 1, 
with the following average properties: burning 
time 10 s, effect height 3 m, net explosive content 
(NEC) 32 g. 

The local particle number concentrations were 
measured over minimum 30 minutes at three 
defined measuring points. The experiments were 
carried out for each measuring point separately 
with identical fountains, since only one measuring 
device (Scanning Mobility Particle Sizer, SMPS) 
was available. The measuring points in the 
examination room are shown in Figure 2. The exact 
same points were also used as validation references 
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(monitor points) for the CFD simulations described 
as follows. 


Numerical modeling 


The overall goal of the simulations with ANSYS 
CFX V 14.0 was to calculate the liberation of smoke 
due to the burn-off of a theatrical pyrotechnic 
fountain and the following transport and dispersion 
of that smoke in the examination room. 


Governing transport equations 


The following equations for the conservation of 
mass, momentum and energy generally apply 
for three-dimensional, transient, compressible, 
laminar Newton fluid flows (in Cartesian tensor 
notation): 


Continuity equation (conservation of mass): 


ap , (pu) i 
Ot Ox, 


Momentum equation: Navier-Stokes (conservation 
of momentum): 


Ol pu, OT. 
OBtb (p J_.@, “1 + pg, 
ot Ox, Ox; OX, 
and 

Ou, Ou; ) 2 du, 
_=n| —-+— |-—7—6, 
ime a cues) al 


Species transport 


OPM, 4 OpUu, LM, = OS. a 
Ot Ox, On. 


and 
Ol, 
J, = PD, als 
Ox, 
and 
m 
Hi =K : 
Ym, 
k=1 


Heat equation (conservation of energy): 


K 
Oph " Opu;h _ 0 A OT) oO Ih, 
Ot Ox,  Ox,\ Ox, ) OX, ta 
ro ye Claes 
Ot Oy, OR, 
and 
K 
h=> 1h, h, =I, + | ¢,,4T 
k=l and (T) 


Ox, 
Nomentlatire 


Latin letters 


Greek letters 


Cp: Specific heat capacity 


g: gravity 

h: specific enthalpy 

i: tensor index (1, 2, 3) 
J: tensor index (1, 2, 3) 
k: species 

m,: Mass 

p: pressure 

S;,: volumic heat flow 
t: time 

T: temperature 

u: velocity 


6, Kronecker—Delta 
n: dynamic viscosity 
A: thermal conductivity 
/: mass fraction 

p: density 

7: Shear stress tensor 


x: Cartesian coordinate 
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Turbulence closure model 


Within the numerical simulations of this work, a 
transient turbulent flow was considered by applying 
the method of Scale-Adaptive-Simulation (SAS). 
This approach uses a dynamic function for the 
turbulent length scale to switch between Reynolds- 
averaged Navier-Stokes equations (RANS) and 
the more sophisticated Large-Eddy-Simulation 
(LES). In contrast to RANS, where all turbulent 
scales are modeled, LES resolves the large scales 
(“eddies”) and models only small scales. The 
final resolution of the turbulence is influenced by 
the mesh size and the time step. SAS basically 
combines the advantages of both approaches: 
reasonable CPU times (in case of RANS) and 


1.500 


3.000 


an increased turbulence resolution and accuracy 
(in case of LES). Further details are given in the 
ANSYS CFX-Solver Modeling Guide.° 


Geometry and numerical grid 


The examination room which was used in 
the investigations of Dutschke ef al.’ was 
implemented into the ANSYS-workbench as a 
“x t’-file (parasolid), generated with a common 
computer-aided design (CAD) tool. The red 
sphere in Figure 2 indicates the location of the 
source point simulating the theatrical pyrotechnic 
fountain. The relevant ventilation openings are 
located near the ceiling and close to the ground. 
Eleven outflow ventilation openings are located in 
the ceiling. Below each bench a row of circular 


TW \\\SYs 
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Figure 3. General view of the domain with two cut planes visualizing the mesh. 
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inflow ventilation openings are located. And, thus 
difficult to visualize in Figure 2 because of its 
small size, two door clearances with a height of 
1 cm and a width of | m are respectively placed on 
the upper floor and the ground floor intersection 
with the wall on the back right side (in the positive 
y-direction) when looking at Figure 2. 


The entire geometry was internally meshed with 
tetrahedrons (unstructured) and prisms for an 
appropriate wallboundary layerresolution, resulting 
in a total number of approximately 3.7 x 10° nodes 
and 12.77 x 10° elements (consisting of 8.73 x 10° 
tetrahedrons and 4.04 x 10° prisms). 


In Figure 3 a general view of the fluid domain is 
illustrated with two planes AA and BB showing the 
mesh. In Figure 4 a detailed view of the AA and 
BB planes details the refinement in relevant areas. 
The mesh was refined around the source point (up 
to the maximum effect height of 3 m), since there 
were high velocity, temperature and concentration 
gradients expected, as well as around the inlet and 
outlet areas. In the same way refinements were 
made in areas with smaller geometry length scales, 
as for example around the tables and benches. 
These refinements related to length scale changes 
around the tables and benches and the prism 
layers for the near wall resolution are shown in 
Figure 5. 


Material properties 


Air as an ideal gas was defined as the existing 
ambient material. The density of the air was 


Plane AA JW Ass 


calculated by the ideal gas law, whereas fixed 
values of 28.96 kg kmol ' for the molar mass and 
1004.4 J (kg K) | for the heat capacity at constant 
pressure were used. The smoke has been accounted 
for as an additional variable with a density of 
1.3 kg m °. The advantage of using an additional 
variable is that no other physical properties have to 
be defined, although a dispersion of the substance 
is simulated. 


Simplifying model assumptions 


To keep the simulations as simple as possible 
in a first step, physical aerosol effects such 
as coagulation (particles collide and coalesce 
with each other, resulting in decreasing particle 
number concentrations), condensation of vapor 
phase on particle surfaces (directly affecting the 
aerosol particle size distribution) and nucleation 
were neglected. Further details on the numerical 
simulation of these effects are given amongst 
others by Hussein et al.* 


The body of the fountain consisting of a cylindrical 
shape of around 3 cm diameter and 10 cm height 
was not explicitly resolved by the CAD model. Due 
to the small size of the fountain body compared to 
the dimensions of the room, neglecting that object 
seemed possible without significant effects on the 
flow in the room aligned with the advantage of 
simplifying the mesh generation and saving CPU 
time. Therefore, instead of modeling the body of 
the fountain, the smoke source was considered as 
a source point situated at the height corresponding 
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Figure 4. Cut planes AA and BB for detailed view of the mesh. 
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Figure 5. Detail side view of the mesh around tables and benches with length scale dependent refinement 


and prism layers. 


to the length of the real fountain. 


Furthermore, other minor simplifications of the 
room geometry were performed for the same 
reasons of making the mesh generation easier and 
saving significant CPU time. 


Boundary and initial conditions 


As the boundary conditions during the experiments 
did not vary much and were therefore highly 
reproducible, mean values were used for the 
simulations. Walls, floor, ceiling, tables and 
benches were modeled as “wall” boundaries 
with a “no slip” condition. The source point was 
defined for releasing a mass fraction of smoke = 1 
with a release velocity oriented in the z-direction 
(normal to the ground in direction of the ceiling) 
of 1.5 ms |. The latter value corresponds to the 
measured velocity during the test runs. The mass 
flow of smoke released was calculated from the 
known values of released volume and release 
duration. 8 L of smoke were released during 10 s, 
corresponding to a volumetric flow of 0.8 L s‘'. 
With an assumed smoke density of 1.3 kg m°, a 
mass flow of 0.001 043 kg s"' was determined. The 
release duration was limited to 10 s according to 
the real functioning time of the fountain. The air 
inlet surfaces were modeled as “inlet” boundary 
conditions, with 0.75 m s ' flow velocity at the 
circular inlets on the ground and 1.39 ms! at 
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the door clearances. The air oulet surfaces on the 
ceiling were set as “outlet” boundary conditions 
with an outflow velocity of 0.87 m s!. The 
velocity values were also calculated from known 
ventilation and geometric conditions. 


As the released smoke is initially at a much higher 
temperature than the ambient air, the heat transport 
has been taken into account using the Total Energy 
model of CFX (see also ANSYS CFX-Solver 
Modeling Guide’ for further details). 


The simulation was set up as a transient simulation 
with a total duration of 1860 s. In the first 60 s 
only the flow in the domain was calculated without 
any release of smoke, to obtain a flow field in the 
domain. After that the release was activated for 
10s and during the last 1790 s the dispersion of the 
smoke was calculated. The time steps used were 
manually preset as a time step list, with variable 
time steps varying from 107 s up to 2.5 10's. 


Processor properties and software 


All CFD simulations were carried out on a high- 
performance cluster consisting of several nodes. 
The relevant node contains of 32 cores in total 
(4 x Magny-Core [8 cores]; AMD Opteron 
2.6 GHz/256 GB RAM), of which 24 cores were 
actually used for the calculations of this work. 
The program ANSYS 14.0 with its workbench 
platform was chosen for the CFD simulations of 
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this work. 


The total wall clock time of one case calculation 
(simulated time of 30 minutes) was around 
14 days. 


Results and discussion 


Although the experimental values were determined 
in three separate experiments, the simulated 
values originate from one simulation run. The 
results shown in the following originate from 
a simulation with the above described mesh of 
nearly 13 million cells. Preceding short estimations 
of the grid influence showed that a highly refined 
mesh is needed. Unfortunately, further refinements 
were in principle possible but not on the available 
hardware. As mentioned before, the concentration 
of smoke during the experiments was measured at 
three different locations in the theatre. In Figure 6 
the measured time dependent concentrations of 
smoke are plotted against the simulated values for 
all three locations. 


Although a perfect match of the transient 
concentrations between the experimental values 
and the simulations is not reached, it can be seen 
that the transient dispersion is reproduced quite 
accurately by the simulations. With only small 
deviations in time, the simulations predict the rise 
of the smoke concentration in good agreement 
with the experimental data. In addition to that, 
the simulation results detailed in Figure 6 also 
reveal the decline of the aerosol concentrations 
due to the comparatively short functioning time 
of the fountain in combination with the existing 
ventilation (boundary) conditions. The overall 
trend is in fair agreement with the experimental 
results. 


Table | gives detailed information about the 


maximum aerosol concentrations and _ the 
corresponding time values of the simulations in 
comparison with the experimental results (see also 
Figure 6). 


The simulations generally tend to over-predict the 
aerosol concentrations in point | and 2 by a factor 
of roughly 6 and 9, whereas in point 3 an under- 
prediction by a factor of 5 was observed. However, 
the simulated times to the respective maximum 
aerosol concentration peaks correspond very well 
with the experimental results (Table 1). 


In terms of labor protection aspects the time 
dependent dose information of a certain aerosol 
exposure of the persons involved is a major key 
factor when assessing the potential hazards of 
the burn-off of indoor pyrotechnics. Therefore, 
the simulated time averaged concentration values 
over 30 minutes were compared at all three points 
with the corresponding experimental results 
published by Dutschke et a/.? As shown in Table 2 
the agreement of the time averaged concentration 
values between simulations and experiments for 
the points | and 2 is better than for the single time 
dependent maximum values presented in Table 1. 
However, for point 3 slightly less good agreement 
on that matter was observed. 


The general trend of over-predicting the aerosol 
concentrations in the close range and under- 
estimating the concentrations further away from 
the fountain was not influenced by the time 
averaging. 


The fact that the general prediction type changes 
from a conservative character at close range to 
the pyrotechnic article to an optimistic quality 
further downstream in the examination room 
may be due to the constant summation of 
calculation uncertainties, implicitly included in 


Table 1 Transient aerosol concentration characteristics 


Points Maximum aerosol concentration [mg m *] Time to maximum concentration [s] 
Experiments* Simulations Experiments* Simulations 

1 1.6 9.4 480 436 

2 (1st relevant peak) 0.11 1.2 480 481 

2 (2nd relevant peak) 0.13 0.8 720 700 

3 0.14 0.028 840 870 


* Values taken from Dutschke ef al.* 
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Figure 6. Measured and simulated aerosol concentrations vs. time at all three measuring/monitor 


points. 
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Table 2. Time averaged smoke concentrations 


es (over 30 min) [mg m *]* 
1 0.39 
2 0.07 
3 0.09 


Experimental time averaged concentration 


Simulated time averaged concentration 
(over 30 min) [mg m *] 


1.62 
0.141 
0.0115 


7 Values taken from Dutschke ef al 


the grid refinement, boundary conditions settings 
and other simplifying model assumptions as 
mentioned earlier in the text. Using mean values 
for the boundary conditions instead of the exact 
values during the experiments could have led to 
differences in the flow and concentration fields. 
Furthermore, although the experimental venue 
was a closed room with regulated ventilation 
conditions, the simulation omits the fact that one 
half of the walls of the room consisted of windows. 
These certainly had an influence on the air and 
aerosol flow conditions inside the room due to 
possible effects of sun radiation. In addition to 


Velocit 
Vector 


1.500e+000 


1.125e+000 


that, the windows might act as heat bridges to the 
outside, leading to differences in the temperature 
field influencing the buoyancy and therefore the 
flow field. 


Figure 7 shows the flow vectors during the burn- 
off of the theatrical fountain 5 seconds after 
ignition on the centerline cross-section area of the 
simulated room. 


The flow inside the room strongly depends on 
several aspects like constructional details (e.g. 
arrangement and dimensions of the chairs and 
tables) and ventilation conditions etc. Three- 
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Figure 7. Flow field in the simulated theatre during the burn-off of a theatrical fountain 


5 seconds after ignition. 
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dimensional flow measurements inside the room 
(e.g. with hot-wire anemometry instruments) 
during the experiments might have given precise 
flow information in order to improve the estimation 
of possible calculation uncertainties. 


The expansion of the aerosol cloud (for a 
concentration of 1 mg m°? and higher) due to 
the dispersion inside the room, as well as the 
contraction due to the short functioning time of the 
fountain and the existing ventilation conditions, is 
displayed in Figure 8. 


Summary 


This study presents results of the investigations 
to evaluate the possibility of using common 
CFD tools for the assessment of possible hazards 
aligned with the liberation of smoke during 
the burn-off of pyrotechnic articles in indoor 
venues. Therefore, a CFD model from ANSYS 
CFX was taken to simulate the dispersion of 
aerosols released by the functioning of a typical 
theatrical pyrotechnic fountain. The calculations 
were afterwards compared with corresponding 
preceding experiments. 


As a general evaluation it can be stated that the 
chosen CFD model assumptions led to overall 
satisfying results compared with the experimental 
values. The calculations over-predicted the 
maximum aerosol concentrations in the close range 
to the pyrotechnic fountain, which corresponds to 
a conservative hazard estimation. The close range 
to the aerosol source point plays an important 
role from the safety point of view, since actors 
and audience are often in very close proximity to 
the pyrotechnic effects used in theaters and stage 
shows. However, a different situation was observed 
in the far range. With increasing distance from the 
source point the maximum aerosol concentration 
values were under-estimated by the simulations in 
this case, likely due to the summation of calculation 
uncertainties corresponding to the chosen model 
assumptions. 


In comparison with the accuracy of the maximum 
concentration values, the transient behavior of the 
aerosol liberation was much better reproduced by 
the simulations. The times at which the respective 
maximum concentration values appeared at all 
three measuring/monitor points matched the 
experimental values very well. 
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Figure 8. Time dependent dispersion of the 
aerosol cloud visualized for a concentration 
of 1 mg m? and higher. 
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Future simulations might even get better 
agreements with experimental results, when: 


e they are in addition validated with flow 
measurements, 


* more sophisticated numerical methods and 
approaches are applied, 


e physical aerosol effects such as coagulation, 
condensation of vapor phase on particle 
surfaces and nucleation are taken into account, 
and 


* amore refined grid is used. 


Conclusions 


In spite of the high numerical efforts, combined 
with comparatively long computation times, an 
application of CFD to estimate potential hazards 
during the use of pyrotechnics with regards to 
the liberation of solid reaction products appears 
to be possible and reasonable. Especially when 
considering the fast development of computational 
processing power inthe recent years, the application 
of CFD in this field will become more and more 
feasible in the future. 


CFD results may also play an important role in 
countries where quantified risk assessments 
(QRA) are taken into account prior to the 
permission of the use of pyrotechnics in stage and 
theatre shows, as they can give information not 
only on absolute concentration values, but also 
on transient dispersion effects, basically for every 
desired venue. These maximum or time averaged 
concentration values could be compared with 
national threshold values with regards to labor 
protection regulations. In particular when using 
several pyrotechnic articles in parallel during a 
show, the advantages of CFD may dominate the 
long computation times. 


Another benefit of CFD is the fact that this 
tool offers the possibility of receiving relevant 
information, such as concentration, temperature, 
velocity etc. for every desired location within the 
3D fluid domain. 


However, several disadvantages of CFD 
simulations will still exist in the future, regardless 
of the promising CPU development. These 
include amongst others the problem or challenge 
of adequately modeling the flow domain (e.g. 
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construction and import of CAD files), and setting 
boundary conditions and material properties 
realistically. 


References 


1 ‘Directive 2007/23/EC of the European 
Parliament and of the Council of 23 May 
2007 on the placing on the market of 
pyrotechnic articles, Official Journal of the 
European Union, 14.6.2007. 


2  prEN 16256:2012 Standard series for 
Pyrotechnic articles — Theatrical pyrotechnic 
articles, consisting of five parts, CEN/TC 
212 WG3, 2012. 


3. A. Dutschke, C. Lohrer, S. Seeger and 
L. Kurth, “Towards the reaction products 
of indoor fireworks”, Proceedings of the 
11" International Symposium on Fireworks, 
Puerto Vallarta (México), 2009, pp. 66—79. 


4 ~ T. Hussein, H. Korhonen, E. Herrmann, 
K. Hameri, K. E. J. Lehtinen, and 
M. Kulmala, “Emission rates due to 
indoor activities: Indoor aerosol model 
development, evaluation, and applications“, 
Aerosol Science and Technology, Vol. 39, 
Issue 11, 2005, pp. 1111-1127. 


5 ANSYS CFX-Solver Modeling Guide, 
Release 14.0, November 2011. 


Page 21 


Recycling of Global Consumer Fireworks 


Andrew Tang, Andy Y. F. Tang and Paul Z. J. Peng 
Tian Cheng Pyrotechnics Laboratory, Lihua Village, Yanxi Town, Liuyang City, Hunan, China 410304 
Email: atang@tcpyrolab.org 


Abstract: Consumer fireworks is one of the two major fireworks sectors, professional display fireworks and 
consumer fireworks, in the fireworks industry. It is estimated that consumer fireworks occupy about half of the 
market share in the global fireworks industry. In terms of packaging materials being used in the consumer 
fireworks sector, it is definitely true that far more material is being used than in the professional display fireworks 
sector. This study is a market survey of global consumer fireworks that are on sale in different countries. It was 
carried out to see what material is being employed in global consumer fireworks and how it is used. The 
environmental pollution from shooting consumer fireworks is considered from different points of view. Recycling 
of consumer fireworks has been introduced in China. Innovation in new materials for consumer fireworks is 
beginning in China. 


Keywords: Consumer fireworks, recycling, packaging materials, environmental pollution, breakdown of 
material, NEC, paper 


Introduction 


Recycling of fireworks is rarely heard of by 
consumers or even people working in the fireworks 
industry. It is commonly accepted that all fireworks 
are used once and disposed of immediately 
afterwards, and that it does not matter how and where 
it is disposed of. In China this is very similar to all 
types of packaging for consumer goods such as 
foodstuffs, unless you want to use the packaging 
again. 


The Green Power organization in Hong Kong has 
conducted an annual survey on people’s consumption 
and celebration habits.’ The analysis shows the habits 
of moon cake wastage, and that recycling moon cake 
packages is still a waste problem. Million of moon 
cake packages are dumped after the Mid-Autumn 
Festival. From study data in 2011 September, there 
were about 47 packaging components used in some 
types of retail moon cake while the average is 11 
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Figure 1. Structure of consumer fireworks (e.g. cake) showing general inner construction. 
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Table 1. Distribution of consumer fireworks samples in the study 


Type No of samples Range 
Total mass 
Aerial shells 5 533-606 
Cakes 15 110-10637 
Firecrackers 9 1.6-857 
Fountains 20 16.5—1944 
Ground spinner 4 10.6—-46.5 
Helicopters 9 10.3—281.5 
Magic whip 1 248 
Missile 3 95-388 
Novelties 7 2.6-217 
Party poppers 5 2.6—-354 
Rockets 5 80.3-359 
Roman candles 5) 70-740 
Sparklers 4 1.3-17.7 
Others 7 : 
TOTAL 99 — 


(g) NEC (g) 
116-133 


16-1124 
0.4-107 
0.4-296 
2.1-8.4 
3.3-67.4 
134 
18.5—70.4 
1.1-200 
0.01-0.7 
16.6—-150 
6.6-107 
0.8-6.9 


a 


“ Some samples were specially designed and so could not be grouped into different groups. The range of gross weight and 


NEC were not meaningful. 


packaging components. 

Normally packaging of consumer fireworks is simple 
and straightforward. The first ancient Chinese 
firecrackers contained no packaging, but were simply 
black powder filled in a piece of bamboo’ closed at 
one end. Packaging for the old style fireworks was 
concerned only with safety and attractiveness. When 
one looked at a firing site after the consumer had used 
fireworks, there were torn papers and broken 
components spread around on the ground. Most of the 
environmental pollution complaints arising from 
shooting fireworks concerned air and noise pollution, 
and much less focused on the left over waste paper 
and components. This study does not describe the 
impact caused by pyro-chemicals but focuses on the 
packing components of consumer fireworks or other 
material that is left over on the ground after using 
fireworks. A preliminary similar report was 
presented . 


Consumption of consumer fireworks 
All consumer fireworks function using different 
formulations, combinations and percentages of black 
powder, stars and effect charges. All of these 
formulas are composed of chemicals. The functioning 
of fireworks creates the effects that are its primary 
performance. As is well known the effects are light, 
gases and sounds that are produced in the air.* During 
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functioning, these chemicals react and generate 
oxides? in the form of small particles, gases, moisture 
and some burnt debris. These small particles or gases 
are blown away by the surrounding air. After the 
functioning of fireworks, all that is left in the air are 
smoke and gases. There is increasing attention to the 
environmental impact of these smokes and gases. A 
lot of studies have been carried out to understand the 
impact and much research and development is still 
continuing to eliminate it. Figure | shows a typical 
construction of one type of consumer firework called 
a battery of shot tubes or the generic term, cakes. 


After functioning of fireworks, all that is left over on 
the ground is packing materials that used to hold 
chemicals. Since people are allowed to handle and 
buy consumer fireworks, the main criteria for 
fireworks packaging are the safety and attractiveness 
of the fireworks. In the early days of fireworks 
manufacturing, all packing materials were made 
mainly of paper and clay. The first firecracker was 
made of black powder and bamboo. Nowadays 
plastics and metals have been introduced. Some 
common items of consumer fireworks are shown in 
Figures 9 to 12. 


Waste from consumer fireworks 
Consumer fireworks are always packed or packaged 
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Figure 2. Percentage of different materials used in con- 


sumer fireworks. 


in an attractive manner to draw consumers’ attention. 
Of course manufacturers will ensure the package 
contains the maximum possible protection to 
consumers so as to provide safe transportation, 
handling and functioning. 


The sole functioning of fireworks generated by pyro- 
chemicals or net explosive content (NEC) is always 
the main concern of both consumer and audience. 
The NEC is composed of chemicals, mainly 
oxidizers, fuels, and reducing agents. These 
chemicals are simple and non-toxic or less toxic 
chemicals such as nitrates, perchlorates, metal 
powders, sulphur, and some metallic oxides. They are 
confined inside the paper tube, shell or both to create 
effects based on different requirements. After ignition 
of the leading fuse and the subsequent functioning of 
the firework, all chemicals react to form different 
forms of oxides. The by-products are water, and 
energy in forms of light and sound, and smoke. The 
residue of fireworks that is left on the ground is the 
waste or containers that used to hold those chemicals. 


It is commonly found in China that most used 


fireworks that are left after firing are collected by 
waste packers because the material is worth money, 
similar to other recyclable materials such as glass, 
metal (aluminum) and plastics. Paper is the most 
collectible material in recycling of fireworks. It is 
collected almost immediately after firing in China. 
Then it is sold to recycling facilities and delivered to 
paper milling factories where paper is manufactured 
and supplied back to fireworks factories. See Figures 
4 to 7 showing fireworks waste collected in China. 


This study focused on the left over waste after the 
chemicals (NEC) have reacted to form gases and 
smoke. 


Experimental 

There were 99 samples randomly selected from 
different factories manufacturing consumer fireworks 
for different markets such as the US, EU, China etc. 
Among these there were different types of consumer 
fireworks, see Table 1. Their gross weights were 
measured and they were dismantled. Different 
materials such as chemicals (NEC), paper, clay, 
plastic, wood and metal were separated then their 
weights measured individually. The majority of these 
consumer fireworks were cakes and fountains 
because they were commonly available in factories 
which implies their large market share. Others were 
rockets, roman candles, small aerial shells, sparklers, 
helicopters, magic whips and missiles. Their 
distribution is shown in Table 1. 


Results and discussion 


Overall average of breakdown of materials in 
global consumer fireworks 

The overall average percentage breakdown of the 
materials used in global consumer fireworks was 
44.6% (w/w) of paper material, 20.7% (w/w) of clay 


Table 2. Breakdown of materials in consumer fireworks by regional market 


Region NEC (%) Pa (%) 
European 20.8 42.2 
Max. 63.8 73.5 
Min. 0.1 22.5 
PRC 16.9 49.4 
Max. 32.0 70.3 
Min. 9.8 30.2 
N. America 21.5 44.3 
Max. 92.0 82.9 
Min. 0.1 2.7 


Pl (%) Cl (%) Wo (%) Me (%) 
5.9 27.5 0.4 28 
66.9 70.0 9.6 36.2 
0.1 2.0 9.6 31.6 
a 32.0 = _ 
74 54.0 = _ 
0.1 7.0 = = 
ibe, 16.0 3.1 1.6 
64.3 64.7 36.4 57.9 
0.2 2.0 2.8 2.0 


Note: No max./min. if there is no average result and absence was not counted as minimum. 
NEC = Net Explosive Content; Pa = Paper; Pl = Plastic; Cl = Clay; Wo = Wood; Me = Metal. 
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Table 3. Breakdown of materials in different types of consumer fireworks 


Type NEC (%) Pa (%) P1(%) Cl (%) Wo (%) Me (%) 

Aerial shells 21.6 59.4 4.0 14.7 0.0 0.0 
Cakes 14.3 39.6 1.7 39.5 4.8 0.0 
Firecrackers 14.4 33.1 1.4 50.1 0.0 0.0 
Fountains 225 49.0 6.7 21.1 0.0 0.0 
Ground spinner 27.9 52.4 1.2 17.7 0.0 0.0 
Helicopters 23.3 56.0 4.5 14.8 0.0 0.2 
Magic whip 54.1 21.8 24.2 0.0 0.0 0.0 
Missile 20.1 21.7 56.0 2.1 0.0 0.0 
Novelties 31.0 42.6 13.6 4.6 5.2 2.4 
Party poppers 0.2 46.6 52.7 0.0 0.0 0.0 
Rockets 28.3 44.9 7.1 8.5 10.7 0.0 
Roman candles 13.6 72.4 0.0 12.0 1.8 0.0 
Sparklers 42.6 19.0 1.3 0.0 0.0 36.9 


material, 9.7% (w/w) of plastic material, 2.1% (w/w) 
of wood material and 1.7% (w/w) of metal material. 
Chemicals (NEC) occupied 20.7% by weight. See 
Figure 2. 


Breakdown of materials in consumer fireworks by 
regional areas 
If these samples were grouped by different regional 


areas such as European countries, USA and PRC, the 
distribution of material was as shown in Table 2. The 
percentage NEC was highest, 21.5%, for samples for 
the North American market, while that for the PRC 
market was the lowest, 16.9%. By contrast, the 
percentage of paper was highest for the PRC market, 
49.4%, and lowest for the European market, 42.3%. 
In terms of clay material, the PRC market was 


Table 4. Breakdown of materials in different consumer fireworks (excluding NEC from the calculation) 


Type NEC Pa-LO (%) 
Aerial shells 0 75.8 
Cakes 0 46.2 
Firecrackers 0 38.7 
Fountains 0 63.2 
Ground spinner 0 Te) 
Helicopters 0 73.0 
Magic whip 0 47.5 
Missile 0 27.2 
Novelties 0 61.7 
Party poppers 0 46.7 
Rockets 0 62.6 
Roman candles 0 83.8 
Sparklers 0 33.1 
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PILLO(%) Cl-LO(%) Wo-LO(%) = Me-LO (%) 
5.1 18.8 0.0 0.0 
2.0 46.1 5.6 0.0 
1.6 58.5 0.0 0.0 
8.6 472 0.0 0.0 
17 24.5 0.0 0.0 
5.9 19.3 0.0 0.3 

52.7 0.0 0.0 0.0 
70.1 2.6 0.0 0.0 
19.7 6.7 75 3.5 
52.8 0.0 0.0 0.0 
9.9 11.9 14.9 0.0 
0.0 13.9 2.1 0.0 
23 0.0 0.0 64.3 


Journal of Pyrotechnics, Issue 31 2012 


Table 5. Breakdown of materials used in cakes 


Materials Average, grams Max. Min 
Gross weight 3376 (100%) 10 637 110 
NEC 397.6 (14.3%) 1124 15.6 
Paper 1459 (39.6%) 4949 38.4 
Plastic 7.3 (1.7%) 48.1 5.1 
Clay 1305 (39.5%) 4644 38.4 
Wood 207.1 (4.8%) 1400 99.0 
Metal 0 0 0 

highest, 32.0% and the North American market was Cakes 


lowest, 16.0%. This shows that fireworks sold for the 
PRC market contain high percentages of paper and 
clay but least NEC. It was envisaged that the local 
transportation of fireworks was less expensive than 
overseas transportation. The packaging of consumer 
fireworks for the PRC market seems to be more 
bulky than that for other markets. 


Breakdown of materials of different types of 
consumer fireworks 

Different types of consumer fireworks function 
differently and so they are manufactured using 
different processes. The materials used, therefore, 
contain different amounts of different materials. For 
example sparklers contain a high percentage of metal 
because of the iron stick although the firecracker does 
not use metal at all. Hence the percentage breakdown 
of material is better presented in terms of its type. 
The breakdown of materials of each type of consumer 
fireworks is shown in Table 3. 

After functioning of these fireworks, all pyro- 
chemicals reacted giving energetic effects such as 
sound, light, gases, moisture etc. The leftover 
materials were paper, plastic etc. Therefore 
percentages of leftover materials were different and 
were calculated without the percentage of NEC. The 
results are shown in Table 4. 


Table 6. Breakdown of materials used in fountains 


Materials Average, grams 
Gross weight 375.8 (100%) 
NEC 70.9 (22.5%) 
Paper 171.9 (49.0%) 
Plastic 13.1 (6.7%) 
Clay 115.9 (21.1%) 
Wood 3.3 (0.3%) 
Metal 0.0 (0.0%) 
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This type was the largest group of consumer 
fireworks available in all markets. Due to the 
similarity of their outer appearance, the sample size 
was smaller than that for fountains. The number of 
shot tubes within each sample differed so the 
performance effect was different too. Hence the 
higher the number of shot tubes and the higher its 
effect, the greater the NEC content of the sample. 
Besides NEC, paper and clay occupied about 80% by 
weight of this type. The breakdown of materials 
included mainly paper (39.6%), clay (39.5%) and 
NEC (14.3%). For a detailed breakdown refer to 
Figure 3. The maximum and minimum weights of 
materials are shown in Table 5. 


Fountains 

The second largest type of consumer fireworks was 
fountains, which extended from a single tube up to 
multiple tubes (battery of shot tubes). Due to the 
varied outer appearance, a larger number of samples 
was randomly chosen. Besides NEC, paper and clay 
were still the two main materials used to construct 
this type of consumer fireworks. More explicitly 
paper occupied about 50% by weight of its gross 
weight. The average waste from fountains consisted 
of paper 49.0%, clay 21.1% and plastic material 6.7% 
by weight. For a detailed breakdown refer to Figure 
3. The maximum and minimum weights of material 
are shown in Table 6. 


Max. Min. 
1,944 40.0 
296.3 0.4 
788.7 12.4 
110.7 0.7 
865.4 2.1 
65.2 0 
1.8 0 
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Table 7. Breakdown of materials in helicopters 


Materials, grams Average Max. Min 
Gross weight 116.7 (100%) 282 10.3 
NEC 24.8 (23.3%) 67.4 3.3 
Paper 64.3 (56.0%) 175.5 6.2 
Plastic 4.8 (4.5%) 15.6 0.02 
Clay 21.3 (14.8%) 79.2 0.7 
Wood 0 (0%) 0 0 
Metal 0.6 (0.2%) 6.1 0 
Table 8. Breakdown of materials in novelties 
Gross Weight 114.4 (100%) 217 2.6 
NEC 37.5 (31.0%) 200 1.1 
Paper 58.5 (42.6%) 160 0.5 
Plastic 9.5 (13.6%) 37.9 0 
Clay 5.7 (4.6%) 19.5 0 
Wood 0.1 (5.2%) 0.9 0 
Metal 2.3 (2.4%) 8.6 0 
Table 9. Breakdown of materials in firecrackers 
Materials, grams Average Max. Min. 
Gross weight 109.8 (100%) 857 1.6 
NEC 13.2 (14.4%) 106.6 0.4 
Paper 63.5 (33.1%) 530.4 0.4 
Plastic 2.2 (1.4%) 18.7 0 
Clay 29.8 (50.1%) 191 0.8 
Wood 0 (0%) 0 0 
Metal 0 (0%) 0 0 
Helicopters the type was more attractive to younger consumers on 


This type of consumer fireworks contained different 
designs. Metal and plastic materials were also 
introduced. There were 9 different samples randomly 
selected from factories. Besides NEC, the first two 
materials were paper, 56%, and 14.8% of clay. The 
lowest was plastic, 4.5% by weight. For a detailed 
breakdown refer to Figure 3. The maximum and 
minimum weights of material are shown in Table 7. 


Novelties 


This type contained all 5 different materials because 
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account of their appearance rather than their 
pyrotechnic effects. There were 7 different samples 
randomly selected from factories. Besides NEC, the 
first two main materials were paper, occupying 
42.6% by weight and plastic, 13.6% by weight. For a 
detailed breakdown refer to Figure 3. The maximum 
and minimum weights of material are shown in Table 
8. 


Firecrackers 


This type of consumer fireworks was the traditional 
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Table 10. Breakdown material of rockets 


Materials, grams Average Max. Min 
Gross weight 192.0 (100%) 359 80 
NEC 56.0 (28.3%) 150 16.6 
Paper 84.5 (44.9%) 195.3 34.0 
Plastic 15.1 (7.1%) 58.2 0.6 
Clay 21.1 (8.5%) 84.4 4.1 
Wood 14.4 (10.7%) 34.4 0 
Metal 0 (0%) 0 0 
Table 11. Breakdown material of small aerial shells 

Materials, grams Average Max. Min 
Gross weight 573.1 (100%) 606.5 533.0 
NEC 123.6 (21.6%) 133.4 116.2 
Paper 340.4 (59.4%) 382.4 313.5 
Plastic 22.9 (4.0%) 28.4 20.5 
Clay 84.3 (14.7%) 98.1 64.3 
Wood 0 (0%) 0 0 
Metal 0.2 (0.1%) 0.3 0 


type. Due to the manufacturing process they did not 
vary too much, and the materials used were mostly 
the same as before. There were 9 different samples 
randomly selected from factories. This type of 
consumer fireworks was separated into batteries of 
firecrackers (or firecrackers on strings) and individual 
firecrackers. Therefore there was a large spread in the 
results. Besides NEC, the first two main materials 
were paper, occupying 33.1% by weight and clay, 
50.1% by weight. This group of consumer fireworks 
did not contain wood nor metal materials. However it 
contained the highest percentage of clay among all 
types of fireworks. It was envisaged that firecrackers 
provided simply a noise effect, a bang, which was 
mainly generated by using clay confinement. For a 
detailed breakdown refer to Figure 3. The maximum 
and minimum weights of material are shown in Table 
9. 


Rockets 


There were 5 different samples randomly selected 
from factories. Besides NEC, the first two main 
materials were paper, occupying 44.9% by weight 
and wood, 10.7% by weight. This type contained the 
highest percentage of wood among all types, because 
all the rocket sticks were of wood. For a detailed 
breakdown refer to Figure 3. The maximum and 
minimum weights of material are shown in Table 10. 
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Small aerial shells 


This type of consumer fireworks was small 
reloadable aerial shells of less than 1.75 inches 
diameter. They could be repeatedly loaded into a 
launch tube and fired one after another. The item is 
commonly available in the American market. Besides 
NEC, the first two main materials were paper, 
occupying 59.4% by weight and clay, 14.7% by 
weight. For a detailed breakdown refer to Figure 3. 
The maximum and minimum weights of material are 
shown in Table 11. 


Roman candles 


There were 5 different sample randomly selected 
from factories for this type of consumer fireworks. 
Besides NEC, the main two materials were paper, 
occupying 72.4%, and clay, 12.0% by weight. For a 
detailed breakdown refer to Figure 3. The maximum 
and minimum weights of material are also shown in 
Table 12. 


Party poppers 


This type of consumer fireworks contained 5 different 
samples but were all very similar to each other. 
Besides NEC, the main two materials were plastic, 
occupying 52.7% by weight and paper, 46.6% by 
weight. For a detailed breakdown refer to Figure 3. 
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Table 12. Breakdown of materials in roman candles 


Materials Average, grams Max. Min. 
Gross weight 310.2 (100%) 740 70 
NEC 44.3 (13.6%) 107 6.6 
Paper 226.2 (72.4%) 544 50.0 
Plastic 0 (0%) 0 0 
Clay 37.1 (12.0%) 88.5 2.9 
Wood 2.0 (1.8%) 1.8 0 
Metal 0 (0%) 0 0 


Table 13. Breakdown of materials in party poppers 


Materials Average, grams Max. Min. 
Gross weight 127.3 (100%) 354 2.6 
NEC 0.3 (0.2%) 0.7 0.01 
Paper 61.1 (46.6%) 173 1.0 
Plastic 65.7 (52.7%) 180 1.6 
Clay 0 (0%) 0 0 
Wood 0 (0%) 1.8 0 
Metal 0 (0%) 0 0 


Table 14. Breakdown of materials in sparklers 


Materials Average, grams Max. Min. 
Gross weight 12.7 (100%) 17.7 1.3 
NEC 4.6 (42.6%) 6.6 0.8 
Paper 3.1 (19.0%) 6.9 0 
Plastic 0.2 (1.3%) 0.8 0.1 
Clay 0 (0%) 0 0 
Wood 0 (0%) 0 0 
Metal 4.7 (36.9%) 9.5 0.5 


Table 15. Breakdown of materials in missiles 


Materials Average, grams Max. Min. 
Gross weight 242.7 (100%) 387.6 95.0 
NEC 48.2 (20.1%) 70.4 18.5 
Paper 48.9 (21.7%) 64.2 23.0 
Plastic 139.0 (56.0%) 243.0 53.3 
Clay 6.2 (2.1%) 9.6 0 
Wood 0 (0%) 0 0 
Metal 0 (0%) 0 0 
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The maximum and minimum weights of material are 
also shown in Table 13. 


Sparklers 


This type of consumer fireworks contained 4 
different samples. It was the type of consumer 
fireworks containing the largest amount of metal 
material. Besides NEC, the main two materials were 
metal, occupying 36.9% by weight, and paper, 19.0% 
by weight. For a detailed breakdown refer to Figure 
3. The maximum and minimum weights of material 
are shown in Table 14. 


Missiles 


This type of consumer fireworks contains 3 different 
samples only. Besides NEC, the main two materials 
were plastic, occupying 56.0% by weight and paper, 
21.7% by weight. For a detailed breakdown refer to 
Figure 3. The maximum and minimum weights of 
material are shown in Table 15. 


Magic whip 


This is the only group containing | sample only, 
because such samples are not common. Besides 
NEC, the main two materials were plastic, occupying 
24.2% by weight and paper, 21.8% by weight. A 
detailed breakdown is shown in Figure 3. 


Conclusion 


Most of the waste from consumer fireworks consisted 
of paper and clay materials after functioning which 
were 44.6% and 20.7% respectively by weight. If 
only leftover materials were counted (deducting 
NEC, 20.7%) these two materials occupied 56.2% 
and 26.1%, a total of 82.3%. Wood occupies 6% and 
metal occupies 3%. All of them are recyclable easily. 
Therefore most consumer firework packaging (waste 
from fireworks) is collected by pickers who sell to 
recycling facilities, see Figures 4 to 7, for recycling 
of paper. Plastic is the most arguable material as it 
depends on which type of plastic it is. The percentage 
of plastic material is about 9.7% by weight (or 9% by 
weight of leftovers) among all consumer fireworks. 
Among these recycling materials, the one with the 
greatest impact on environmental pollution is plastic 
material. Most of the plastic found in the study was 
PVC and PE. 


If packaging used more environmentally friendly 
materials, the percentage of recycling would increase, 
and it would have less impact on the environment. 
Some factories have started to develop the use of new 
reusable materials in the manufacturing of fireworks, 
such as plastic tubes for reloadable aerial shells or 
some new materials which are biodegradable in the 
environment, as shown in Figure 8. As a general 
trend, manufacturers have become aware of the 
challenge as well as the opportunity offered by the 
environmental impact caused by fireworks. 
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Figure 3. Material used in different types of consumer fireworks by percentage. 
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Figure 5. Trash paper from different packaging con- Figure 6. Recycling paper by mixing trash paper and 
sumer products. fibers before milling. 
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Figure 7. Finished recycled paper. 


Figure 8. Innovative material to replace paper and 


clay in manufacturing cakes 


Figure 9. Some examples of consumer fireworks: (a) missile; (b) helicopter; (c) firecrackers; (d) novelty. 


Table 16. Breakdown of materials in magic whip 


Materials 
Gross weight 
NEC 

Paper 

Plastic 

Clay 

Wood 

Metal 


Average, grams 
247.5 (100%) 
133.8 (54.1%) 
53.9 (21.8%) 
59.8 (24.2%) 
0(0%) 

0 (0%) 

0 (0%) 


Max. 
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Abstract: Dynamic, non-isothermal analysis of the pyrotechnic mixture of a sound-producing fire 
cracker, KNO;/AIl/S/H3BO3, was carried out under Nz atmosphere by using a_ simultaneous 
thermogravimetric and differential thermal analyser. Intermediates and the residue formed in each stage 
of thermal decomposition of the pyrotechnic mixture were characterised by X-ray diffraction analysis 
and the most probable reaction pathway was proposed. Thermogravimetric (TG) and differential 
thermogravimetric data (DTA) were used for the interpretation of the mechanisms and kinetics of 
decomposition by means of a model-fitting method, the Coats—Redfern equation, and model-free 
methods, the Arrhenius and Kissinger methods. The values of activation energy (E) were calculated as 
11-16 kJ mol', 24-53 kJmol' and 15-56kJ mol' for stages I, II and III respectively by the model- 
fitting method, with the Avrami—Erofe’ev model (A2), two-dimensional model, nucleation and growth 
mechanism. The pre-exponential factor (In A) of each stage of thermal decomposition at various linear 
heating rates was calculated and the probable decomposition mechanism was proposed for the first 
decomposition step as the expulsion of sulphur, and the second and third decomposition steps as the 
conversion of aluminium to alumina are obtained. Similarly, the reaction was also found to be first 
order. 


Keywords: pyrotechnics, fireworks, thermal analysis, activation energy, Coats—Redfern equation. 


directly as thermal energy or more usually as 
light, sound and kinetic energy to achieve a 
desired effect. This reaction must not begin to 
proceed as soon as the pyrotechnic composition 
is mixed and the spontaneity is prevented by the 
characteristic property of pyrotechnic 
compositions, the ‘activation energy’ barrier. An 
externally applied flame can initiate the ignition 
process by supplying the activation energy to the 
composition. If the activation energy barrier is 


Introduction 


Pyrotechnic compositions are mixtures of 
reducing and oxidising chemicals that are 
capable of undergoing self-sustained 
combustion. The reducing agent is often referred 
to as fuel. In general, the composition of 
fireworks is a mixture of sulphur or phosphorus, 
perchlorate or nitrate and pure aluminium 
powder.’ The pyrotechnic mixtures are sensitive 


to impact, friction, electrostatic thermal energy 
and auto-ignition temperature. The composition 
of fireworks should possess high sensitiveness to 
flame and low sensitiveness to other factors like 
impact, friction and electrostatic energy. 


A pyrochemical reaction  characteristically 
produces heat energy which may be useful 


low, accidental ignition will be more likely 
because a small amount of mechanical or 
electrostatic energy can cause ignition of the 
composition. Different oxidizers are used in the 
pyrotechnic mixture in order to get the desired 
effect of the fireworks and the parameter, energy 
of activation, can be used to check which 
oxidizer in the pyrotechnic mixture can be 
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effective for a particular type of fireworks. 


Potassium permanganate and phosphorus have 
low activation energy and therefore could not be 
used in the composition of pyrotechnics in the 
firework industry. Activation energy is not only 
affected by the nature of the oxidizers and the 
igniter but the same materials in different 
combinations may also lead to low activation 
energy creating a significant hazard. Accidental 
ignition during processing or spontaneous 
ignition during storage is the major problem in 
the fireworks industry which is the main focus of 
our work. The mechanism involved in the 
combustion and _ explosion’ reaction of 
pyrotechnic mixture in a firework is entirely 
different compared to the spontaneous ignition 
reaction of pyrotechnic mixture during storage. 


In the firework industry, after mixing the 
ingredients, the pyrotechnic mixture is used to 
fill contrivances. Cardboard boxes are used for 
cake bombs and atom bombs; tubes are used for 
Chinese crackers and cones are used for flower 
pots and the combustion reaction is carried out in 
a confined way in the fire crackers. Differential 
Scanning Calorimetry (DSC) and Accelerating 
Rate Calorimetry (ARC) are suitable techniques 
for calculating thermodynamic parameters in a 
closed system as well as to interpret the 
mechanism. Simultaneous thermal analysis (TA), 
thermogravimetric (TG) and differential thermal 
analysis (DTA) can be used for studying the 
decomposition kinetics of the pyrotechnic 
mixture for ignition or explosion during storage. 


The parameters that influence the thermal 
behaviour are particle size distribution, purity 
and moisture content of the chemicals etc. 
Decrease in particle size of KNO; had an adverse 
effect on sensitivity to both impact and friction to 
a greater extent than the particle size of Al,” and 
set the lowest onset temperature for the 
exothermic activity. The heat content from the 
exothermic activity was also found to vary with 
particle size.” Decreasing the particle size from 
micro-size to nano-size pyrotechnic mixture was 
found to produce sound effectively and it was 
reported that a smaller amount of pyrotechnic 
mixture is required to produce the allowed sound 
level of 120- 125 dB(A) as prescribed by the 
Govt. of India* on using a nano pyrotechnic 
mixture.’ Although it is expected that the particle 
sizes create large interfacial areas and increased 
number of atoms at the particle interface which 
on ignition lead to a higher heat of reaction, the 
various physical and chemical processes 
occurring concomitantly and competitively may 
affect the process. 
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The effect of particle size of the pyrotechnic 
mixture on the heat of reaction was carried out 
using DSC° and ARC.’ Many studies have been 
carried out on the factors affecting the sound 
level produced from the , fire crackers’ and 
especially on particle size** which played a key 
role in producing the sound level from the fire 
crackers. Kinetic parameters for all the stages of 
decomposition have not been reported so far. In 
the present work, thermal analysis was carried 
out to deduce the decomposition pattern of the 
pyrotechnic mixture under well controlled 
conditions by TGA and DTA techniques. Kinetic 
parameters for different decomposition stages are 
derived from the results and the conclusions 
drawn from them are presented. Importance is 
given to calculating the kinetic parameters, 
energy of activation (£) and pre-exponential 
factor (In A) by following model-fitting methods, 
the Coats-Redfern equation, and model-free 
methods, the Arrhenius and Kissinger methods. 
Data for energy of activation related to different 
proportions of pyrotechnic composition and for 
commonly used mixtures of oxidizers and igniter 
have not been reported in the literature. An 
attempt is made to calculate the energy of 
activation for pyrotechnic mixtures with different 
oxidizers in varying compositions. 


The present work focuses on the most commonly 
used oxidizer in sound producing fire crackers, 
KNO; with sulphur, aluminium and boric acid. 


Experimental 
Chemicals and materials 


The chemicals used for the preparation of fire 
cracker were obtained from a_ firework 
manufacturing company situated in Tamilnadu, 
India. The purity and assay of the chemicals are: 
KNO; — 97.6%, Sulphur (S) — 99.9%, 
Aluminium (Al) — 99.8% and boric acid (H3BO3) 
— 99%. Aluminium powders were of grade 999 
(200 mesh, 75 microns), KNO; of 120 mesh and 
125 microns, S of 100 mesh and 150 microns and 
H3BO3 of 100 mesh and 150 micron sizes were 
used. All these chemicals were sieved through a 
100 mesh brass sieve. The samples were stored 
away from light and moisture. All these samples 
were mixed in the ratio of KNO;3: Al: S : H3BO3 
as 57.5/20/22/0.5 %. 


Instruments 
Thermal analyser 


Thermal analysis (TA), thermogravimetric (TG) 
and differential thermal analysis (DTA) were 
carried out using a Perkin-Elmer, Pyris diamond 
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model thermal analyser with a heating rate of 1, 
10 and 20°C min", and a temperature range of 
the standard system of room temperature to 
900 °C. 


X-Ray diffraction analyser 


The pyrotechnic mixture and its intermediate 
products on decomposition were analysed by X- 
ray analysis by employing a Siemens 800 X-ray 
diffractometer (XRD) with Cu Ka radiation with 
a view to characterising the intermediate 
products of decomposition. Thermal runs of the 
pyrotechnic mixture were interrupted in the 
TGA-DTA system (heating rate 10K min’, 
cooling rate 60 K min ') and the intermediates in 
the range of temperatures, 280-300 °C and 950- 
980 °C were collected and the collected residues 
were analysed by XRD. Adequate care was 
taken to avoid moisture pick up before XRD 
analysis. 


Results and discussion 


Figure 1 shows the simultaneous TGA, DTG and 
DTA curves of pyrotechnic mixture obtained at 
a heating rate of 10 Kmin! under N> 
atmosphere. The decomposition patterns 
obtained under other heating rates are similar 
and the increase in heating rate leads to an 
increase in the characteristic peak temperature 
(7,) which can be seen in the thermograms 


DTA, uv 


(Figure 2). 


DTG analysis showed that decomposition 
occurred in three stages: the first stage ranges 
from 150 to 250°C, the second stage is in the 
range from 580 to 680°C and the third stage 
ranges from 750 to 850°C (Table 1). 


Decomposition occurred only in three stages by 
increasing the heating rate. In TG analysis, the 
second and third stages are not clearly 
differentiated. | Pyrotechnic mixture before 
decomposition and residues after the first and 
third stages of decomposition were collected and 
were analysed by XRD analysis (Figure 3). 
Powder diffraction patterns are typically plotted 
as the intensity of the diffracted X-rays versus 
the angle 20. Peaks will appear in the diffraction 
pattern at 26 values when _ constructive 
interference is at a maximum. By measuring the 
20 values for each diffraction peak, the d- 
spacing (the distance between the diffracting 
planes) can be calculated automatically by the 
data analysis software for all of the peaks in the 
diffraction pattern. The 26 values obtained from 
the XRD pattern were compared with standard 
JCPDS (Joint Committee on Powder Diffraction 
Standards) data and marked in the XRD pattern 
in Figures 3-5. 


The sharp mass loss at 150—250°C is assigned to 
the removal of sulphur making a weight loss of 
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Figure 1. 7G, DTA and DTG curves of KNO;/S/Al/H3BO; in N> atmosphere; heating rate: 10K min 1. 
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Figure 2. 7G analysis of KNO;/S/Al/H3BO; in N> atmosphere; at three different B: 1, 10, 20 K. 


$ Pyrotechnic mixture 
1000 # ALS 
$ KNO,* 
S@ 
H,BO, # 
800 * 373 


600 


arb. units 


400 


200 


20 40 60 80 
20 degrees 
Figure 3. XRD pattern of pyrotechnic mixture, KNO;/S/Al/H3;BO3;, before decomposition. 
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Figure 4. XRD pattern of the residue after decomposition between 280 and 300 °C. 
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Figure 5. XRD pattern of the residue after decomposition between 950 and 980 °C. 
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Table 1 Comparison of the expected and the observed mass loss for the decomposition of KNO3/Al/S/ 


H;BO; at 10 K min‘ under nitrogen atmosphere 


Weight loss (TG data) (%) 


Stages Temperature Assumed reactions 
Observed Expected 

S —S()—- S(g) 

I 150-250 °C 24 22 
2KNO; — 2KNO) + O2(g) 

I 580-680 °C 15 15.84 
4Al + 30,(g) 2 ALO; 

Ul 750-850 °C 35 35.29 


4KNO, — 2K,0 + 4NO(g) + O2(g) 


24% by exothermic reaction. The exothermic 
reaction indicated the conversion of sulphur into 
molten melt by melting or fusion and is a 
physical process that results in the phase 
transition of a substance from a solid to a liquid, 
at which the ordering of ionic or molecular 
entities in the solid breaks down to a less ordered 
state and the solid liquefies. Substances in the 
molten state generally have reduced viscosity 
with elevated temperature; an exception is the 
element sulphur, whose viscosity increases to a 
point due to polymerization and then decreases 
with higher temperatures in its molten state. 
Sulphur is not removed by sublimation which is 
an endothermic phase transition of a substance 
directly from the solid phase to the gas phase 
without passing through an intermediate phase 
and requires additional energy. 


In the fire crackers, it is assumed that all the 
sulphur taken in the pyrotechnic mixture is 
converted to sulphur dioxide and the reaction is 
highly exothermic providing energy to initiate 
the reaction of the pyrotechnic mixture. The 
absence of peaks due to sulphur (Figure 4) 
confirms the removal of sulphur in the first stage 
as shown in the XRD of the residue at stage I. 


The second stage of mass loss which occurred 
between 580-680°C is _ assigned to 
decomposition of oxidizer, KNO3;, liberating 
oxygen. Liberation of oxygen from the oxidizer 
causes 15% weight loss and the reaction is 
exothermic. The third stage of decomposition 
occurred immediately after stage II of 
decomposition between 750-850°C with 35% 
mass loss corresponding to the expulsion of NO 
gas and as the heating rate increases, the 
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temperature difference between these two stages 
is narrowed down. Complete exothermic reaction 
of pyrotechnic mixture occurred at stage III 
leaving only Al,O3 and KO as residue which 
were collected for XRD analysis. The XRD 
pattern (Figure 5) supported the suggested 
mechanism and confirmed the final residue 
obtained from the fire cracker is only potasstum 
oxide and aluminium oxide (Table 1). 


Based on the observations of weight loss from 
TG data of pyrotechnic mixture by three stages, 
the mechanism given in Table | is assumed to 
take place. There is good agreement between the 
expected theoretical mass loss and the observed 
mass loss for all the stages of decomposition of 
pyrotechnic mixture of KNO;: Al: S: H3BO3 in 
the ratio of 57.5/20/22/0.5 % shown in Table 1. 


The kinetics of the thermal decomposition of 
sound producing pyrotechnic mixture, KNO:/S/ 
Al/H3BO3 under N> atmosphere at different 
heating rates, were studied by model-dependent 
and model-free methods for each _ stage 
separately. 


Determination of energy of activation, E,: model- 
free methods 


The transformation rate during a reaction is the 
product of two functions, one depending solely 
on the temperature, 7, and the other depending 
solely on the fraction transformed, a: 


da/dt =f (a).K(T) (1) 


where da/dt is the derivative of the fraction 
converted with respect to time. It was calculated 
for every 10% mass loss of the mixture for all 
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Figure 6. Fraction reacted, a, versus temperature for stages I, II and III at three different heating 


rates. 


the three stages of decomposition in the 
temperature range of 150—250°C, 580-680°C 
and 750-850°C respectively at different heating 
rates. a 1s defined by the expression as 


a=(%w i- %w,)/(Yow; — Yows) (2) 


where %v; is the mass percent at any time ¢ and 
%w;, and Yow, are the initial and final mass 
percent sample,” ’° respectively. A typical plot of 
a versus temperature for stages I, II and III is 
shown in Figure 6. 


The temperature dependent function is generally 
assumed to follow an Arrhenius type’ ” 
dependency: 


Ink = InA—E,/RT (3) 


Thus, to describe the progress of the reaction for 
non-isothermal experiments at all temperatures 
and for all temperature-time programmes, the 
reaction rate at all times depends on both f(a) 
and A(T), and hence the determination of f(a), 
In A and E (the so-called kinetic triplet) are 
needed. The non-isothermal DTG data (Figure 7) 
were used to calculate the kinetic parameters 
using Arrhenius method. The linear Arrhenius 
plots of Ink versus 1/T for the first stage of 
decomposition of pyrotechnic mixture at various 
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heating rates were plotted and similar type plots 
were obtained for other stages also (Figure 8). 


From the slope, the activation energy (E,)'° for 
the decomposition of the pyrotechnic mixture 
was calculated (Table 2). The activation energy 
values obtained are 30.7 + 3.4, 58.23 + 4.2, 
58.83 + 4.6 kJ mol! for stages I, II and III 
respectively. 


The activation energy for the non-isothermal 
decomposition of pyrotechnic mixture was also 
calculated from the TG data using Kissinger 
expression: 


In (B/Tm’) = In [n(1 — dm)" 'AR/E,] — E,/RTm (4) 


where fis the heating rate, A is the pre- 
exponential factor, F, is the energy of activation, 
T,, and G@, are the absolute temperature and mass 
loss at the maximum mass loss rate (da/df),, and 
R is the gas constant. This method yielded values 
of 29.9 + 2.0, 45.1 + 2.0 and 108.1 + 2.0 kJ mol! 
for stages I, II and III, respectively (Table 2) 
from the slope of In (6/ ie ) as a function of 1/7; 
at the maximum mass loss rate (Figure 9). 


The fraction reacted, a, 0.10—0.9, was also used 
for the kinetic analysis using the Flynn—Wall 
method'® at different heating rates £ using the 
expression: 
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Figure 7. DTG analysis of KNO;/S/Al/H3BO; in N> atmosphere at three different heating rates. 
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Figure 8. Arrhenius plot of Ink versus temperature for non-isothermal decomposition of KNO;/S/Al/ 
H3BO; in N2 atmosphere; heating rate: 1,10 and 20 K min *. 
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Figure 9. Kissinger plot for stages I, II and III in N atmosphere for the non-isothermal decomposition 
of KNO;/S/Al/H;BO; using TG data. 
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Figure 10. Apparent activation energy versus a, fraction reacted, for stages I, II and III in Nz atmosphere 
for non-isothermal decomposition of KNO;/S/Al/H3BO3  . 
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Figure 11. Mass loss (%) for stages I, IT and III versus temperature for the decomposition of KNO;/S/ 
Al/H;BO; at different p. 
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Figure 12. 42-mechanistic plot of In g(a)]T° versus 1000/T for KNO;/Al/S/H3BO; in N> atmosphere at 
different P. 
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Figure 13. FJ plot of [—In(1 — a)] versus 1000/T by KNO;/Al/S/H3BO; in Nz atmosphere 
at different P. 


Table 2 Activation energy for the decomposition of KNO;/Al/S/H3;BO;3 


Methods Stages E/kJ mot ' Correlation coefficient (7) 
I 29.945 —0.96452 
Kissinger II 45.145 —0.97819 
Il 106.0+2 —0.98159 
I 30.744 —0.89253 
Arrhentus I 58.2+6 0.92733 
Il 58.93 +6 —0.9157 
I 10.644 —0.95209 
iso-conversional i 17.046 0.82797 
Ill 60+ 10 —0.92403 
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Table 3 List of solid-state rate equations used in the present study 


Description 


Rate controlling mechanism 

P1 power law 

E1 exponential law 

A2 Avrami—Erofe’ev equation | . 
A3 Avrami—Erofe’ev equation 2 
A4 Avrami—Erofe’ev equation 3 
B1 Prout-Tompkins 

R2 contracting area 

R3 contracting volume 

D1 one-dimensional diffusion 
D2 two-dimensional diffusion 


D3 three-dimensional diffusion, Jander’s 


D4 three-dimensional diffusion, Ginstling Brounshtein 


F1 first order 
F2 second order 


F3 third order 


In B = In (AB,/R) — In[F(a)]—E,/RT (5) 


The activation energy for each conversion point 
(E,, a) was calculated from the slope of the linear 
plot of the Flynn—Wall method at different 
heating rates f. The plots of E, versus a, fraction 
reacted for stage I, II, III (Figure 10) show that 
E, decreases with the extent of conversion (a), 
0.2 <a < 0.9 and the activation energy is high at 
the initial fraction reacted which indicates that 
the decomposed products are in equilibrium with 
the gaseous reactants. The decreasing of FE, with 
increasing extent of conversion indicates that 
porous type solids decomposed more easily. The 
results reveal that the dependence of apparent 
activation energy (F,) on the extent of 
conversion helps not only to disclose the 
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Equation 
g(a) 

lt 

Ina 

[-In(1— a)]"” 
[-ind=a)]”* 
[-In(1- a)]"" 


In[a/(1— a)] 


1-(1-a)'” 
1-(1-a)'? 
oe 


(l-a)ln(l-a)+a 
M-(d-a'3P 


(1 — 2a/3) - (1 - a)? 


—In(1 - a) 
1/1 — a) 
[V/ — a)? 


complexity of a decomposition process, but to 
identify its kinetic scheme as well. 


The E, value is lower for stage I and II (10.444 
and 15.88 kJ mol ') compared to the Z, value of 
the third stage of decomposition 
(50.25 kJ mol '). 


Dissociation mechanism — model-fitting method 


Kinetic analysis by model-fitting of the thermal 
data obtained for all the three stages of 
decomposition was done in order to arrive at the 
mechanism of the dissociation processes and to 
deduce the kinetic parameters. The rate of the 
reaction under non-isothermal conditions can be 
expressed by the following relation: 
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Table 4 Kinetic parameters and possible rate controlling processes of decomposition of KNO;/AI/S/ 


H3BO3 
TG data 
Stage Mechanism BK min! see: 
E/kJ mol! In A/min coethicient' (7) 
Stage I A2 1 16.054 3.662 0.96027 
10 12.771 ~4.950 ~0.97311 
20 11.916 ~6.778 ~0.97636 
Stage Il A2 I 53.316 38.02 ~0.88901 
10 25.884 4.518 ~0.92033 
20 24.195 2.423 ~0.90682 
Stage III A2 1 56.222 3.253 ~0.94685 
10 15.729 ~7.516 ~0.7897 
20 17.212 5.868 ~0.90674 
da/dT = K(T)/b — fla) (6) line when the correct g(a) function is used in the 


where a is the fraction reacted at temperature 7. 
fia) is the conversion function which is 
dependent on the mechanism of the reaction. / is 
the rate of heating employed in the experiment 
and k(T) is the rate constant as a function of 
temperature. Equation (6) can be represented by 
its integral form as follows: 


In[g(a)/T"] = In[(AR/BE)(1 -2RTIE)] — E/RT (7) 


The algebraic expression of the integral g(a) 
functions that are tested in this work are listed in 
Table 3. A plot of In g(a)/T° versus 1/T gives a 
straight line when the correct g(a) function is 
used in the equation. The g(a) function describes 
the mechanism of the reaction. Straight lines 
with a high correlation coefficient and low 
standard deviation were selected to represent the 
possible controlling mechanism. The plot of 
mass loss for the stages I, II and III with respect 
to the temperature is given in Figure 11. The 
fraction reacted, a, was evaluated as a function 
of temperature. 


The kinetics of all the three stages of the thermal 
decomposition of sound producing pyrotechnic 
mixture, KNO;/S/Al/H3BO; in N> atmosphere at 
different heating rates was followed by 
employing the Coats—-Redfern’’ approximation 
(Table 3) which gives the expression. 


A plot of Ing(a)/T° versus 1/T gives a straight 


Journal of Pyrotechnics, Issue 31 2012 


equation. The g(a) function describes the 
mechanism of the reaction. Straight lines with 
high correlation coefficient and low standard 
deviation were selected to represent the possible 
controlling mechanism. The corresponding 
kinetic parameters were then calculated and are 
shown in Table 4. The model-fitting method is 
performed in the conversion region where the 
apparent activation energy is approximately 
constant and where a single model may fit. 


The non-isothermal kinetic data of mixture in the 
limit of 0.2 < a < 0.9 are fitted to each of the 
thirteen reaction models which are indicated in 
Table 3. The values of activation energy (£,), pre 
-exponential factor (In A) and the coefficient of 
linear correlation (7) for various kinetic models 
at three different heating rates are presented in 
Table 4. The Arrhenius parameters (£,, In A) are 
highly variable, exhibiting a strong dependence 
on the reaction model chosen. On the other hand, 
more than one model, namely A2 and A3, gave 
fairly good coefficients of linear equations. 
However by comparing the FE, value with the 
model-free method, the A2 mechanism 
(Figure 11) was accepted. 


The corresponding kinetic parameters were then 
calculated and are shown in Table 4. The best fit 
for the first decomposition step of expulsion of 
sulphur, and the second and third decomposition 
steps of converting aluminium to alumina are 
obtained using A2, the two-dimensional Avrami 
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Table 5 Kinetic parameters and order for KNO3/Al/S/H3BO; 


TG data F 

Stage Order BIK min™ ener india cocffictent (2) 
Stage I Fl 1 24.72276 21.60864 0.96452 

10 18.25368 4.58985 ~0.97819 

20 16.55692 0.9595 ~0.98159 
Stage II Fl 1 100.05091 91.57722 0.89253 

10 45.28849 24.68004 ~0.92733 

20 41.92852 20.51068 ~0.9157 
Stage III Fl 1 46.16195 22.35522 —0.95209 

10 25.29485 0.93511 ~0.82797 

20 28.22128 4.19184 -0.92403 
—Erofe’ev model which corresponds to kinetically complex. The order of the reactions 


nucleation and growth mechanism. Similarly, the 
order of the reaction was also found to obey F1 
(Figure 12), a first order reaction and the 
calculated kinetic parameters are given in 
Table 5. 


Conclusions 


Thermal analysis of pyrotechnic mixture, KNO3;/ 
Al/S/H3BOs, from sound producing fire crackers 
was carried out at three different heating rates, 1, 
10 and 20K min’. The thermal decomposition 
of the pyrotechnic mixture occurred in three 
stages. The pathway of the decomposition was 
predicted by analysing the XRD pattern of the 
residue at stages I and III as sulphur in the 
mixture is first decomposed exothermically, and 
the liberated heat energy is used to initiate 
further exothermic reactions in the mixture 
leaving K,O and Al,O; as final residue. The 
kinetics of thermal decomposition were studied 
for all three stages and the energy of activation, 
E,, was calculated as 10-30 kJ mol ', 20-58 kJ 
mol! and 56-70 kJ mol ' for stages I, II and III 
respectively by the model-free methods, 
Arrhenius, Kissinger and _ isoconversional 
methods which are in good agreement with the 
energy of activation, EF,, calculated as 11- 
16 kJ mol ', 24-53 kJ mol! and 15-56 kJ mol ! 
for stages I, II and III respectively by the model- 
fitting method A2, a two-dimensional model, 
nucleation and growth mechanism. The 
dependence of effective activation energy on 
extent of conversion shows that the process is 
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obeys first order with the energy of activation, 
E,, as 16-24 kJ mol ', 41-53 kJ mol! and 25- 
46 kJ mol ' for stages I, II and III respectively. It 
is our duty to identity some chemicals which can 
be added to the pyrotechnic mixture so as to 
increase the E, value to prevent autoignition of 
pyrotechnic mixtures while in storage in the 
fireworks industry. The calculation of the 
kinetic factor E, will help to identify the best 
composition of pyrotechnic mixture for storage. 
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